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Abstract 

Ten conservation laws in useful polynomial form are derived from a Cartan form and Exterior 
Differential System (EDS) for the tetrad equations of vacuum general relativity. The Noether 
construction of conservation laws for well posed EDSs is introduced first, and an illustration given, 
deriving 15 conservation laws of the free field Maxwell equations from symmetries of its EDS. The 
Maxwell EDS and tetrad gravity EDS have parallel structures, with their numbers of dependent 
variables, numbers of generating 2-forms and generating 3-forms, and Cartan character tables all 
in the ratio of 1 to 4. They have 10 corresponding symmetries with the same Lorentz algebra, and 
10 corresponding conservation laws. 
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I. INTRODUCTION 



In a recent paper [1| on the mathematical structure of the WEBB tetrad partial dif- 
ferential equations of vacuum relativity, they were expressed as an Exterior Differential 
System (EDS) in a space of 44 dimensions. These equations had been formulated for nu- 
merical integration and the EDS was set with variables labeled to conform with tetrad 
notation, as 4 to-be-independent coordinates x l , 16 fields that become, in any solution, 
orthonormal tetrad components l Xj, and 24 Ricci Rotation coefficients Tijk(= —Tn~j), but 
all treated equally ab initio. Guided by Cartan's method of the movable frame, the EDS 
was given in terms of just certain combinations of basis 1-forms, viz. 9i = r]ij J dx k and 
Uij = Tkji k Xi dx l . i,j,= 1...4. rjij is a set of constants expressing the orthonormality rela- 
tions and signature of the tetrad basis in a solution. The EDS was generated by four 2-forms 
(torsions, T % = d6 l + u l j A #■?) and four 3-forms for Ricci-flatness, and, for closure, the four 
additional 3-forms that are the exterior derivatives of the 2-forms. A Cartan character anal- 
ysis of the EDS was reported, demonstrating that it encoded a well-posed system of partial 
differential equations having the Cauchy properties of uniqueness, and evolution from initial 
conditions. Ten invariance generators of the EDS were noted. 

A Cartan 4-form A was found for this EDS giving, when pulled back to a solution, a 
Lagrangian density L for deriving the partial differential equations as the Euler-Lagrange 
equations of a functional variational principle. A Cartan 4-form is unique only up to the 
exterior derivative of an arbitrary 3- form (so-called boundary terms), but with one choice 
the Cartan 4-form found, evaluated on a solution, was the the Ricci scalar or Hilbert action 
density. With a different choice, it evaluated to Det\ k \i\ times a quadratic function of the 
Tijk which we now know is found in various first order functional Lagrangian formalisms 
beginning with early work of Weyl |3| and M0ller 0. Functional variation of L respecting 
the relations (from the EDS) for the "intensity" fields T^k as curls of the "potential" fields 
h \i lead back to the tetrad partial differential equations. 

Weyl was concerned principally with relativistic fields and unification, and there followed 
several long traditions of tetrad and spinor field theories with quadratic Lagrangians for 
generalized gravitation such as Einstein-Cartan theory and Poincare gauge theory; cf. Pel- 
legrini and Plebanski [f| , Hehl [?J , Tupper and Phillips , Hayashi and Nakano t 9] , Hayashi 
and Shirafuji Nester and Tung jll|, Meyer 1^1 ; for considerable subsequent literature 
cf. Itin 0. 

M0ller in Ref. [5] saw tetrad fields as more basic than metric fields g^ v within classical 
Einstein general relativity. From his Lagrangian, he derived an associated energy-momentum 
complex. The M0ller energy-momentum complex for general relativity was clearly discussed 
by Cattaneo at GR4 in 1965 by du Plessis and most recently, with many more cita- 
tions, in a very useful review by Szabados ^6|. These authors emphasize that the complex 
associated to the M0ller Lagrangian is distinct from those of Einstein/Freud, Rosen, Lan- 
dau/Lifschitz, Bergmann/Komar, Goldberg and others (including earlier work by M0ller.) 

In the pure vacuum case of general relativity considered in [1], the novel use of a Car- 
tan form and EDS shows that we in fact have a classic field theory, and its ten invariance 
generators (or gauge transformations) enable us to go on and directly write ten polynomial 
conservation laws for the field equations. Four of the invariances can be understood as com- 
ing from (or describing the possibility of) arbitrary coordinate variation in a solution, and 
the conservation laws to which they lead, when pulled back to a solution, give a polynomial 
expression of the M0ller complex. The other six invariances similarly correspond to arbitrary 
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Lorentz transformation of the tetrad frame at each point of a solution, and their associated 
conservation laws appear to be the (distributed, nonlocal) spin of the vacuum gravity field; 
they all should be useful in validating the integrations of numerical relativity 

In Section II we summarize the EDS version of Noether theory, relating invariances of 
a Cartan form A to conservation laws. In Section III we show as an example how this 
works with an EDS and Cartan form for Maxwell theory, resulting in 15 conservation laws. 
We point out the parallel structure to this of the tetrad EDS, Cartan form and invariance 
generators of Ref. and in Sections IV and V use the latter to obtain the ten polynomial 
conservation laws. 



II. NOETHER THEORY FOR EDS 

The generators of continuous invariance transformations of a well posed EDS have been 
called its isovectors, and a systematic method for their derivation illustrated by examples, 
including source-free electromagnet ism ^?}- An isovector is a vector field V, in the space 
where an EDS is given, generated by a closed exterior differential ideal / (dl C /), that 
satisfies 



Isovectors of EDSs systematize the Lie theory of invariances of systems of partial differential 
equations. We note that an EDS is not unique: sometimes some variables and generating 
forms can together be dropped, or it can be prolonged, in the sense of Cartan, by introduc- 
ing more (jet) variables, equivalent to higher partial derivatives, together with additional 
generating forms in the ideal /. Or, by partially integrating, new variables and generating 
forms can be introduced, potentials and so-called pseudopotentials in prolongation structures 
[18], now often called integrable extensions 0]. The well posed nature of an alternate EDS 
should be checked by calculation of its Cartan character table, and involutory variables. 
The group of isovectors will change, and in particular more isovectors can arise, under such 
prolongation to new but related EDS. 

A conservation law for an EDS is an n-l-form \1/ (n being the number of independent 
variables in a solution, 4 in the EDSs treated below) not in the I generating the EDS, 
but whose exterior derivative dty is in /, and hence vanishes on solutions. By the Stokes' 
theorem the integral of \I/ over a closed 3-boundary in a solution vanishes. Physically, this 
is conservation of a nonlocal quantity that is the integral of \l/ over a given open three 
dimensional domain, for example over a 3- volume taken while holding t = constant. 

Use of a Cartan form, together with knowledge of its associated invariances (that may well 
be isovectors), leads to an elegant version of the Noether theorems for finding conservation 
laws from Lagrangians. We recall the defining property of a Cartan form A for an EDS 



This says that the multisymplectic, or Poincare-Cartan form dA, when it is multiplied by- 
contracted with- any arbitrary vector, say X, yields a 4-form that is in /: 



£ V I c / 



(1) 



dA C I A I 



(2) 



X\dA c / 



(3) 



In the mathematical literature 2(| a Noether vector A" for dA is defined to satisfy 



£ N dA = 



(4) 
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For a particular choice of boundary term for A a somewhat different but natural definition 
of Noether vector could be 

£ N A C I (5) 

Both of these definitions lead to conservation laws. A third, more restrictive, definition of 
Noether vector, which we shall use in the following, satisfies both of the above: 

£ N A = (6) 

The conservation law associated to Noether vector N is simply \I/ = VJA; if \[> is not in / 
expanding equations (5) or (6), together with equation (3), is the required relation 

C / (7) 



III. CONSERVATION LAWS FOR SOURCE-FREE ELECTROMAGNETISM 

In Ref. [l[ it may have erroneously been implied that all previously known multisym- 
plectic forms for Maxwell's equations were sums of terms, each of which was the product 
of a contact 1-form and a 4- form. One cited paper, by Robert Hermann |2l|], gave a new 
multisymplectic form for electromagnetism which was the product of a 2-form and a 3-form. 
The resulting Cartan form and EDS, and conservation laws, provide a simple example of 
Noether theory that is a close parallel to those which arise for the WEBB system to be 
treated in Sections IV and V. For tetrad relativity the numbers of potential and intensity 
variables, numbers of generating 2-forms and 3-forms, number of terms in dA and the Cartan 
characters of successive ranks are all just four times larger than those for this Maxwell EDS. 

To be very explicit, the Maxwell EDS is set in a 14 dimensional space, with variables 
which we shall successively denote (x 1 , x 2 , x 3 , x 4 , Ai, A 2 , A 3 , At, F 23 , Fzi, F\ 2 , -F14, F 24 , F 34 ). 
This order is chosen consistently to allow cycling 1 — > 2 — > 3 — > 1. The EDS is generated 
by a 2-form, its closure 3-form and a second 3-form: 

dAi A dx 1 + dA 2 A dx 2 + dA 3 A dx 3 + dA 4 A dx 4 - 
Fudx 1 A dx 2 + F^dx 1 A dx 3 - F^dx 1 A dx 4 - F 23 dx 2 A dx 3 - F 24 dx 2 A dx 4 - F u dx 3 A dx\8) 

dF u A dx 1 A dx 2 + dF 1A A dx 1 A dx A + dF 23 A dx 2 A dx 3 + 

dF 2 4 A dx 2 A dx 4 - dF 31 A dx 1 A dx 3 + dF 34 A dx 3 A dx 4 (9) 

dF 12 A dx 3 A dx 4 - dF 14 A dx 2 A dx 3 + dF 23 A dx 1 A dx 4 + 

dF 24 A dx 1 A dx 3 + dF 31 A dx 2 A dx 4 - dF 34 A dx 1 A dx 2 (10) 

The Cartan character table of this EDS, in the notation of Ref. HI, is 14(0,1,3,5)4+1; it 
is well posed and the x l are involutory with s 4 = 1 gauge freedom. A 4-form for radiation 
gauge could be added to achieve 14(0,1,3,6)4. (With gauge similarly specialized, the WEBB 
characters are 44(0,4,12,24)4 0.) 

If the 2-form is dropped the two exact 3-forms by themselves generate a simpler well 
posed EDS in 10 variables for just the Maxwell field equations without potentials. The 
character table then is 10(0,0,2,4)4. It has a 17 parameter isogroup [17|. It does not 
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however have a Cartan form. Its simplest prolongation appears to be introduction of 4 
potential fields, giving the EDS we treat here, and allowing Hermann's Cartan form. Two 
of the isovectors are lost (duality rotation and field scale change) but the 15 remaining 
become Noether vectors having the Lie algebra of the conformal group C(3,l). They were 
first found by Harry Bateman in 1908 [i^. We give them explicitly below. The associated 
15 conservation laws were derived from them by Bessel-Hagen j23[, using Noether theory. A 
plethora of further, probably less interesting, symmetries and conservation laws (beginning 
with Lipkin's zilch |24|) follows from other prolongations, to jet variables and successively 
higher partials of the fields. We can only refer to a sample of this literature, to recent papers 
with comprehensive bibliographies |25| [26| [2?} . 

Returning to the EDS in 14 variables, the exterior product of the two 3-forms vanishes 
identically, so it is clear that there is a multisymplectic form, the Hermann form, the outer 
product of the 2-form Equation (8) and 3-form Equation (10) 

dA = dA x A dF 12 A dx 1 A dx 3 A dx 4 - dA 1 A dF u A dx 1 A dx 2 A dx 3 + dA x A dF 31 A dx 1 A dx 2 A dx 4 + 
gL4 2 A dF 12 A dx 2 A dx 3 A dx 4 - dA 2 A dF 23 A dx 1 A dx 2 A dx 4 - dA 2 A dF 24 A dx 1 A dx 2 A dx 3 - 
dA 3 A dF 23 A dx 1 A dx 3 A dx 4 - dA 3 A dF 3± A dx 2 A dx 3 A dx 4 - dA 3 A dF 34 A dx 1 A dx 2 A dx 3 - 
dA 4 A dFu A dx 2 A dx 3 A dx 4 + dA 4 A dF 24 A dx 1 A dx 3 A dx 4 - dA 4 A dF 34 A dx 1 A dx 2 A dx 4 + 

F 12 dF l2 A dx 1 A dx 2 A cfe 3 A dx 4 - F u dF 14 A dx 1 A dx 2 A dx 3 A dx 4 + 
F 23 dF 23 A dx 1 A cfe 2 A dx 3 A dx 4 - F 24 dF 24 A dx 1 A c/x 2 A cfe 3 A dx 4 + 
F 31 dF 31 A dx 1 A c/x 2 A cfe 3 A dx 4 - F 34 dF 34 A dx 1 A dx 2 A cfe 3 A do{il) 

Any integral of this gives a Cartan form. Using boundary terms to eliminate bases such 
as dF\ 2 makes it easier later to pull it back into a solution to write a first order functional 
Lagrangian, so we take 

A = FudAx A dx 1 A dx 2 A dx 3 - F 3 idA x A dx 1 A dx 2 A dx 4 - F 12 dA x A dx 1 A dx 3 A cb 4 + 
F 24 dA 2 A dx 1 A dx 2 A dx 3 + F 23 dA 2 A dx 1 A dx 2 A dx 4 - F 12 dA 2 A cfe 2 A dx 3 A cb 4 + 
F 34 dA 3 A cfe 1 A cfe 2 A dx 3 + F 23 dA 3 A dx 1 A dx 3 A cfe 4 + F 3 idA 3 A <ix 2 A dx 3 A cfe 4 + 
F 34 cM 4 A dx 1 A dx 2 A dx 4 - F 24 dA 4 A cfe 1 A dx 3 A dx 4 + F 14 dA 4 A dx 2 A <ix 3 A <ix 4 + 

2 ( F 2 2 s + ^li + ^i 2 2 - F u - F L ~ F u) dxl A dx2 A rfa;3 A d&2) 

Pulling back this A to a solution manifold, using Equation (8), yields the functional 
Lagrangian as usually given in terms of the fields Ai (actually missing from this Lagrangian) 
and their partials (here in the combinations Ajj — Aj t i). 

The Bessel-Hagen conservation laws of vacuum electromagnetism follow directly from 
contraction of the vectors of the isogroup on this A. For conciseness we will tabulate only 
seven typical isovectors. The other eight can be obtained from these by a simple cycling 
of the coordinate indices 1^2^3^1in both the component labels and entries. The 
vectors denoted L a , L 4 , a = 1, 2, 3 are space-time translations, and lead to well-known energy- 
momentum density expressions for the vacuum electromagnetic field. The P a 's and Q a 's 
satisfy the Lorentz Lie algebra, and give spin and cospin field densities behaving like angular 
momentum and center of mass motion conservation. R a , R 4 , Rq generate the extension of this 
to the conformal group C(3.1), and no intuitive understanding of the associated conservation 
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laws seems to have been offered. 



Li = (f.O, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0) 

L 4 = (0,0,0,1,0,0,0,0,0,0,0,0,0,0) 

P 1 = {x\0,0,x\-A 4 ,0,0,-A u F 2A ,-F u ,0,0,F l2 ,-F 31 ) 

Q, = {0,-x 3 ,x 2 ,0,0,-A 3 ,A 2 ,0,F 31 ,-F l2 ,0,0,-F 3A ,F 24 

R = (x\ x 2 , x 3 , x\ -A t , -A 2 , -A 3 , -At, -2F 12 , -2F 31 , -2F U , -2F 23 , -2F 24 , -2F 34 ) 

i?i = (^(a 42 - a; 2 ' - ^ + x 42 ), x x x 2 , x 1 ^ 3 , x x x 4 , -A x x l - A 2 x 2 - A 3 x 3 - A4X 4 , 

A x x 2 - A 2 x x , A x x z - A3X 1 , -A4X 1 - A ± x 4 , -2F X2 x x + F 23 x 3 + F 24 x 4 , 
-2F3XX 1 + F 23 x 2 - F 34 x 4 , -2FUX 1 - F 24 x 2 - F 34 x 3 , -2F 23 x x - F 31 x 2 - F 12 x 3 , 
-2F24X 1 + F 14 x 2 + F 12 x 4 , -2F34X 1 + F 14 x 3 - F 31 x 4 ) 

-R4 — — ^ X X ? X X , X X ; ~( X X X X ) ; A4.X ~\~ A\X j A4.X ~\~ J\-2X j 

A4X 3 + A3X 4 , A x x x + A 2 x 2 + A 3 x 3 + A4X 4 , -F 24 x x + F 14 x 2 + 2F i2 x 4 , 
F 3i x l - F 14 x 3 + 2F 31 x 4 , F 12 a; 2 - F 31 a; 3 + 2F 14 x 4 , -F 34 x 2 + F 24 x 3 + 2F 23 x 4 , 
-F l2 x l + F 23 x 3 + 2F 24 x 4 , F 31 x l - F 23 x 2 + 2F 34 x 4 ) (13) 



IV. ENERGY AND MOMENTUM CONSERVATION IN TETRAD GRAVITY 

We now consider the A of vacuum tetrad gravity, Equation (14) of Ref. [1]: 

A = uj is A u' 5 A 6 k A e t e ijkl - 2o; y A u ks A 6 s A 0,e iiW + 2u^- A T fc A 0^' fc ' (14) 

A has 1440 terms when written out. The coordinates x l do not explicitly appear in A so it 
is immediate that £q/q x iK = 0, i.e. the Lj are both isovectors and Noether vectors: 

L x = (1,0,0,0,0,....) 
L 2 = (0,1,0,0,0,....) 
L 3 = (0,0,1,0,0,....) 

L 4 = (0,0,0,1,0,....) (15) 

(43 zero entries in each component list!) We have used Mathematica to calculate the four 
conserved 3-forms if. 

U = Li\A (16) 

and gone on to pull them back to a solution (using the curl relations Equation (17) of Ref. 
0]), to there expand them as polynomials on the coordinate 3-form basis: 

t p = T l p dx j A dx k A dx l e ijk i (17) 

The are evidently a non-local energy-momentum complex for gravitation, depending of 
course on the choice of this particular A, and its boundary terms, as has been explained. T % v 
has 1248 terms in each of the 12 off-diagonal components, and 1584 terms in each of the 4 
diagonal components! They are homogeneous polynomials, each term being quartic in the 
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l \j and quadratic in the r*- fc , and can readily be evaluated point-by-point for integration, 
given a numerical solution for these 40 fields. The apparent complexity of can however 
be greatly reduced by using summations over the tetrad indices and by introducing the 
transpose inverse fields i\ j (cf. Ref. Q. These arrays are the dual contravariant tetrad 
basis vectors in solutions.) Perhaps the most succinct way of writing the result is then the 
following: 

T* = Det\%\ i\ q q j X p (18) 

where 

C] = Z) + y) + T-C j +W j (19) 

and 

z) = (i/2)(r^-r fc }r^) 
y) = v kH v lkJ - v%vf 

rtj — i J- jki - i m l j 

n) = (i/2)(r fc M r^-r^r-)5i (20) 

Equations (18-20) are the first principal result of this paper, and express in explicit, cal- 
culable, form the nonlocal conservation laws for energy/momentum found by M0ller (if. 
M0ller also sought to justify adjoining six additional equations to achieve locality; these we 
now recognize as the "Lorentz" gauge conditions of van Putten and Eardley [ij. Without 
restricting the gauge, in Section V we next find six additional nonlocal conservation laws . 



V. SIX SPIN CONSERVATION LAWS FOR TETRAD GRAVITY 

Solutions of the tetrad equations go into neighboring solutions under arbitrary orthogonal 
or Lorentz transformation of the tetrads, that is, keeping the r]ij unchanged. While physically 
identified, these are mathematically distinct, so generated by isovectors of the EDS and 
leading to nontrivial conservation laws. Both A and the generators of the EDS, Equation 
(6) of Ref. [H, are covariant with respect to Lorentz indices of the l Xj and F^k so the 
components of these isovectors can be written by inspection. We will denote them P a and 
Q a as they satisfy the same Lorentz Lie algebra as the Maxwell invariances of Section III. 
First, recall the order of variables adopted in Ref. [l|: 

„.l „2 3 _4 1 \ h 1\ 1\ 2x 2 x 2 x 2 x 3 x 3 x 3 x 3 x 4 x 4 x 4 x 4 x 
X ,X ,X ,X , Ai, A2, A3, A4, Ai, A2, A3, A4, Ai, A2, A3, A4, Al, A2, A3, A4, 

Tll2, rxi3, TlU, Tl23? 1^124, 1^134, 1^212) 1^213) 1^214, 1^223) 1^224, 1^234, 1^312) 1^313) ^14, 

1^323; 1^324, T334, T412, T413? T414, T423; T424, T434 (21) 

The components of two isovectors then express, respectively an infinitesimal tetrad boost in 
the iA* direction 

P x = (0, 0, 0, 0, - 4 A X , - 4 A 2 , - 4 A 3 , - 4 A 4 , 0, 0, 0, 0, 0, 0, 0, 0, - X A 2 , 

1^412 — 1^124; 1^413 ~~ ^134: "414, T423, 

T424 — rn2, T434 _ Tll^i — T224, — ^34, , 0, — P 2 12, — P213, "1^324, 

— r 334 , 0, 0, — r3i2, — L313, rn2 — r424> Fu$ — r 434 , 

FlU, T123; Tl24 — T412, Tl34 — ^13) (22) 
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and rotation about 



Qi 




and the other four isovectors follow from cyclic permutations 1 — > 2 — ► 3 — > 1. We have in 
the last two equations specialized to the tetrad EDS with rjij = diag(+, +, +, — ), which is 
of most interest, and to give correspondence with the Maxwell example. 

Finally it is straightforward to calculate six conserved 3-forms by contracting these vectors 
on A. We label them M pq = —M qp . Each of these conserved 3-forms (without further pull 
back) is already in terms of the four bases dx l A dx^ A dx k , each of whose coefficient is 
a polynomial with 48 terms, third order in the l Xj and linear in the Y^. The result is, 
as in Section IV, most conveniently written using the contravariant fields jA- 5 , and we use 
raised/lowered Lorentz indices to allow for general rjif 



As remarked by Szabados [16], these expressions were first found by M0ller as "superpo- 
tentials" for the tetrad energy /momentum complex, Eq. (18). Like the energy /momentum 
conservation laws, these spin conservation laws should be of use in connection with numerical 
integrations. 
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